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Abstract 

Effective Boltzmannians in the sense of the block spin renormalization 
group are computed for the 2D Ising model. The blocking is done with 
majority and Kadanoff rules for blocks of size 2 by 2. Transfer matrix 
techniques allow the determination of the effective Boltzmannians as 
polynomials in u = exp(4/3) for lattices of up to 4 by 4 blocks. The 
zeroes of these polynomials are computed for all non-equivalent block 
spin configurations. Their distribution in the complex (3 plane reflects 
the regularity structure of the block spin transformation. In the case 
of the Kadanoff rule spurious zeroes approach the positive real (3 axis 
at large values of (5. They might be related to the renormalization 
group pathologies discussed in the literature. 



1 Introduction 



Regularity is at the heart of position space (block spin) renormalization 
group. It is usually assumed, and of central importance, that coupling con- 
stants of the block spin effective Hamiltonian depend in a non-singular way 
on the parameters of the original theory. There are, however, situations 
where this assumption is not valid. Pathological behaviour in renormaliza- 
tion groups of the low temperature Ising model was first observed by Israel [0J 
and Griffith and Pearce || |3| . An extensive and rigorous analysis of regular- 
ity properties and pathologies in Ising model block spin transformations was 
performed by van Enter et al. M. The central observation is that in certain 
situations the effective measure for the block spin theory cannot be repre- 
sented as exp(— H). This means that the effective measure is non-Gibbsian. 
See also for a careful analysis of the situation. 

In this paper, I present some numerical results on the distribution of 
zeroes in the complex j3 plane for block spin Boltzmannians of the 2D Ising 
model. These Boltzmannians are partition functions with "fixed" block spins 
fi, viz. 

B(M)=£P(A*,ff)exphMT(*)]- (1) 

cr) encodes the blocking rule. Why should these zeroes provide interest- 
ing information? If the usual renormalization group assumptions are true, 
the zeroes of B(fj,) should - for all block spin configurations fi - behave "bet- 
ter" than those of the full partition function. Note that the zeroes of the 
full partition function approach the real axis at the critical point ||. This 
should not happen for the zeroes of B(fi)\ Furthermore, one might expect 
that the pathologies described in the literature are related to the distribution 
of zeroes close to the j3 axis at large positive values. 

This article is organized as follows: In section 2 the model notation is set 
up, and the blocking rules are defined. Section 3 gives a sketch of the transfer 
matrix technique used to compute the polynomials. Section 4 summarizes 
the numerical results for the majority rule blocking. Observations for the 
Kadanoff blocking rule are reported in section 5. Conclusions follow. 
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2 Model and Block Spin Definition 



We deal with the 2D Ising model, with partition function 

Z = $>xp[-/3#(<7)] , (2) 

(7 

where 

H(a) = - £ a % o 3 . (3) 

<i,j> 

The CTj assume values ±1 and are defined on a square lattice of extension L 
by L, supplied with periodic boundary conditions. The energy if is a sum 
over all pairs of nearest neighbours. In the infinite volume limit the model 
undergoes a second order phase transition at (3 C = |ln(-y/2 + 1) = 0.4406868. 
A block spin transformation with scale factor 2 is defined as follows. For L 
even, the lattice is divided in blocks of size 2 by 2. Given the configuration of 
the cr-spins in a block f, a block spin fij is chosen with probability P(fi,a). 
The majority rule is defined through 

p {v,v) = II Pi(Pi> ff ), ( 4 ) 

blocks / 



with 



Pl (^j,a) = { 1 if \ii J2 ie i <*i > , (5) 
else . 
The so-called Kadanoff rule is 

exp^/i/Eie/O-i) 

a) = - — — . (6) 

2 cosh (u Eig/O-j) 

In the limit u — > oo one recovers the majority rule. 

Given that Eu -P^) °") = lj the full partition function can be rewritten 

as 

^ = EBM, (7) 

with 

B(//)=£P(//,a)exp[-/?#(<T)]. (8) 

(7 

Usually, one aims at a parametrisation B(fi) = exp[— H' (//)], where if' is the 
effective Hamiltonian. Note, however, that this is impossible in the patho- 
logical situations discussed in the literature Q. 
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3 Computation of B(fi) 



An exact computation of B(fj) as function of j3 seems impossible. However, 
on lattices up to at least L = 8 it can be determined by numerical transfer 
matrix calculations. L = 8 corresponds to 4 by 4 blocks, with 2 16 = 65, 536 
block spin configurations \l. We do not have to consider them all. Con- 
figurations connected through global spin flip or geometric symmetries (re- 
flections, shifts, rotations) have the same effective Boltzmannian B(fi). A 
careful counting yields for the number of non-equivalent configurations N(L) 



4 6 



10 



N(L) 4 13 479 86056 

The result for iV(10) was taken from ref. |7|]. 

After multiplication with a constant pre-factor, B(n) can be expressed 
as a polynomial of order L 2 in u = exp(4/3). The coefficients of this polyno- 
mial can be computed by transfer matrix multiplication. In order to avoid 
notational complication, I give a sketch of the method for the case of 3 by 
3 blocks (L — 6). The generalization to other values of L is then obvious. 
Some of the notation is depicted in table [I]. The effective Boltzmannian can 
be expressed as 

B(n) = Tr [T ■ S(fa,fa,fa) ■ T ■ S(fa, fa, fa) ■ T • S(fa, fa, fa)) . (9) 

Here, S is a 2 6 by 2 6 matrix, labeled by the Ising row configurations. It 
depends explicitly on the line configuration of prescribed block spins. E.g., 
the matrix elements of S(fa, fa, fa) are 



S(fa, fa, fa)a,r = exp 



H Pl (fa,a,r). (10) 



L i=l J 1=1 

The matrix T, which is also of size 2 6 by 2 6 , is defined by 

6 



exp 



L i=l 



(11) 



It is not difficult to represent the transfer matrix multiplications in eq. (|9D 
in terms of operations on the coefficients of polynomials in u. The computer 
implementation of these operations form the basis for the results presented 
in this paper. 
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Table 1: Notation for transfer matrix in the case L = 6. 

4 Results for the Majority Rule 

Let us start with some results for L = 4. The four non-equivalent configura- 
tions, called c#l . . . c#4, are specified in the head of table §. In the columns 
we quote the coefficients Bk(n) of the polynomial B(/jl) = Sfc=o-^fc(/-0 uk - 

The zeroes of these polynomials were determined with the help of the com- 
puter algebra program MapleV. For zeroes uq not lying on the negative real 
axis, we then computed the corresponding (3o values through (3q = |ln(-u )- 
The distributions of these numbers for L = 4 are shown in figure 0, with 
different symbol code for the four block spin configurations. The figure also 
contains a circle of radius f3 c around the origin. Note that the two zeroes 
closest to the critical point belong to configuration c#4. 

The results for the 3 by 3 block lattice are plotted in figure |^. The 
13 non-equivalent block spin configurations are specified in table ^. One 
observes that again the zeroes closest to the critical point belong to the fully 
magnetized configuration (c#13). 

The zeroes of the 479 effective Boltzmannians on the 4 by 4 block lattice 
are shown in figure || The plot also shows (with crosses) the zeroes of the full 
partition function. It seems that the Boltzmannian zeroes do not approach 
the real axis in the critical region, whereas the full partition function zeroes 
do. To check this in more detail, we compare the distribution of zeroes in the 
critical region for the three available lattice sizes together, see figure f|. The 
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Table 2: Coefficients in the polynomials B(p) = Y^k=o Bk(n) u k , for 
the four independent 2 by 2 block spin configurations on an L = 4 
lattice, majority rule. 



5 



0.8 
0.6 
0.4 
0.2 


-0.2 
-0.4 
-0.6 
-0.8 




"c#l" 
"c#2" 
"c#3" 
"c#4" 



-0.6 -0.4 -0.2 0.2 0.4 0.6 0.: 
Re beta 



Figure 1: Zeroes of B(fi) for the four non-equivalent block spin 
configurations on a 2 by 2 block lattice, majority rule. The circle has 
radius [3 C . 
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Figure 2: Zeroes of B(n) for the 13 non-equivalent block spin con- 
figurations on a 3 by 3 block lattice, majority rule. 
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Table 3: The 13 non-equivalent block spin configurations on a 3 by 
3 lattice. 
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Figure 3: Zeroes of B(fi) for the 479 non-equivalent block spin 
configurations on a 4 by 4 block lattice, majority rule. In addition, 
the zeroes of the full partition function on an 8 by 8 lattice are shown 
(crosses) . 



plot clearly demonstrates that the zeroes do not move towards the real axis 
in the critical region. One might conclude from this plot that there should 
exist in the L — > oo limit a strip around the real (3 axis ranging from (3 = 
at least up to (3 C where the effective Boltzmannian is free of (3 zeroes. In 
this region it should thus be possible to take the logarithm without danger. 
Furthermore, high temperature (small (3) expansions for the renormalization 
group could by analytical continuation be used in the critical region. 

There is another observation when comparing figures || and |3|. With 
increasing L, more and more zeroes populate the part of the plane with 
larger real part of (3. They are not obviously approaching the real axis there, 
but we also cannot exclude such a scenario. Note that in the analysis of 
van Enter et al. M the case of the 2D Ising majority rule was not treated. 
It is therefore presently not clear whether in this case a large (3 pathology 
exists. We shall see in the next section that in case of the Kadanoff rule 
(where pathologies do exist) spurious zeroes seem to approach the axis at 




§ : Jm. 
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Figure 4: Common plot of the Boltzmann zeroes for L = 4, 6, and 
8, majority rule. 



large positive 0. 



5 Results for Kadanoff Rule 

It has been observed many times that passing from a ^-function block rule to 
a "Gaussian smeared" rule improves the locality and analyticity properties 
of the effective Hamiltonian, see e.g. ||. In this section some results will 
be presented on the distribution of zeroes in the case of the Kadanoff rule 
eq. (§) for uj = 1 and uj = 2. For finite uj there is a finite probability that the 
block spin does not have the same sign as the majority of spins in the block: 

Sigj^J probp = 1)^=1 prob(/i = 1)^=2 
4 0.99966465 0.99999989 

2 0.98201379 0.99966465 

In figures [5| and |6] we show the zeroes for L = 4 and L = 6, both for 
uj — 1 and oj — 2. Obviously, most of the zeroes in the neighbourhood of the 
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Figure 5: Zeroes of B(fi) for L = 4, for the majority rule (w = oo), 
and for the Kadanoff rules with to = 1 and u = 2. 



critical circle do not move very much. In fact they are already very close to 
their majority rule values. However, compared with the u = oo case, extra 
zeroes appear that populate the region of larger real part of (3. It might be 
interesting to note that the zeroes with the largest real part of (3 belong to 
block spin configuration c^2, followed by those of c#l. Another observation, 
also clearly seen in figure |6], is that the spurious zeroes move to the right when 
uj is increased. Most likely they are shifted to infinity when passing to the 
majority rule. 

In figure 7 the zeroes of all three lattice sizes, L = 4, 6, and 8 are plotted 
for the case w = l.[| A careful inspection reveals that the zeroes in the right 
half plane do move towards the real axis when the lattice size is increased. 
They might thus very well be reflecting the large (3 pathologies. 

^or L = 8, some 16 of the 479 block spin configurations were not taken into account 
because the search for the zeroes of the corresponding polynomials suffered from numerical 
instabilities. 
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Figure 6: Zeroes of B{fj) for L = 6, with Kadanoff rule. Top: uj = 
bottom: u = 2. 
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Figure 7: Zeroes of B(fi) for L = 4, 6, and 8, with Kadanoff rule, 
u = 1. 
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Conclusions 



The distribution of zeroes of effective Boltzmannians was studied for 2D Ising 
systems. Both in the case of the majority and the Kadanoff rule a finite 
region around the critical point stays free of zeroes, also when the volume is 
increased. In case of the Kadanoff rule, however, zeroes populate the right 
half plane and approach the real axis at large (3. They might be related 
to the pathologies of a number of Ising renormalization groups discussed in 
the literature. It would be very interesting to understand this relation (if it 
exists). Furthermore, one should try to understand the origin of the extra 
zeroes that appear for finite u. Van Enter et al. were not able to extend 
their analysis of pathologies to uj = oo in two dimensions. The present study 
shows that the distribution of large (5 zeroes changes significantly when u> 
becomes finite. The present findings do not, however, exclude the possibility 
of large (3 pathologies of the 2D Ising majority rule renormalization group. 
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